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Abstract 

We propose a realization of quantum computing using polarized photons. The 
information is coded in two polarization directions of the photons and two- 
qubit operations are done using conditional Faraday effect. We investigate 
the performance of the system as a computing device. 
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After the early discussion of quantum computing the field has attracted much 

attention because Shor ||^ has shown that the famous factorization problem can, in principle, 
be speeded up considerably by quantum data manipulation techniques. The recent work on 
quantum computations has been reviewed by Bennett in Ref. Many realizations have 
been suggested, at present the most promising seem to be ions trapped electrodynamically 
P or in a cavity p. 

In the recent work 0, one of us considered the possible use of photon polarization 
states to carry quantum information. The advantage is that they provide a natural two- 
state basis with no additional Hilbert space components, such as the vacuum state, that 
may constitute losses of the coding. The single photon coding allows an easy detection, in 
contrast the vacuum state is hard to distinguish from a failed detection. The photon coding 
also allows long dephasing times and the possibility to transfer the information from one 
device to another through fibers. The purpose of this paper is to investigate how realistic 
this suggestion is by numerical integration of a semirealistic situation. 

The elliptically polarized photon state [a^a\ + a_al)|0) can be manipulated by the 
Faraday effect induced by the presence of a second photon h\.. These are supposed to 
selectively transfer population from the ground state |0) in Fig.l to the levels | ±1). Because 
each photon a)^ sees only one transition | ± 1) — |2) it becomes modified by the population 
transferred by the photons 6!^ If we keep the transition | ± 1) ^ |2) off resonance, the atom 
acts as a dielectric only and hence the relative phases of a!^ are modified; this is a turning 
of the axis of the state (a;+a+ + a;_al)|0). It was shown in that this allows the gated 
application of an arbitrary unitary transformation. 

In this paper we are looking at two different cases. Case I corresponds to Fig.l when 
Ai = 0. Here both transitions |0) | ± 1) are in resonance and both photons a\. experince 
a modified phase. The situation is symmetric: if 6+ is present alone we achieve a phase shift 
exactly opposite in sign to that caused by the presence of only. In Case II we detune one 
of the transitions, |0) — | + 1) say; see Fig.l. Then only the presence of the resonant photon 
6+ affects the phase of the a^-photons. This corresponds to the gate where the presence of 
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the state fe-|0) does nothing. Most gates discussed earher in the hterature are of this type. 
The four-level system shown in Fig.l is described by the Hamiltonian 

H = ^2(0^0+ + ala_) + ni{blb+ + blb^) 

+CJ2|2)(2| +CJ1+I + 1)(+1| +cji_| - 1)(-1| 

+a;o|0)(0| + Ai + 1)(0| + 6+| - 1)(0| + h.c.) 

+A2 (a-|2)(-l| + a+|2)(+l| + h.c.) . (1) 

In Case I we assume that the states | ±1) are degenerate and that the transitions |0) | ± 1) 
are at resonance, A{ = uji± — ujq — Qi = 0. The transitions |±1) ^ |2) are assumed detuned, 
i.e. A2 = UJ2 — ujq — Qi — Q2 = — — ^2 is nonzero. In Case II we lift the degeneracy 
of levels I ± 1) by setting uJi^ 7^ uji_. Then the transition |0) ^ | — 1) is taken at resonance 
cui- — luq = Qi but the detunings Al^ = — uq — Qi and = UJ2 — ujq — Qi — Q2 are 
nonzero. The transition | + 1) |2) is detuned by A'2 = LJ2 — uii+ — fl2 = — A(^; this 
is assumed well off resonance too. 

The initial state is taken to be the disentangled form 

|^,„) = (a+at + a-al){f3+bl + (3-bl)\0), (2) 

where |0) denotes the vacuum of the fields. The coefficients are in general complex numbers 
normalized to unity. We propagate the state vector (j^) to the time t with the Hamiltonian 
(|l|) and write the final state as 

|^o.i) = e-^^*|v|/,„)=X:aK), (3) 

1=1 

where we have numbered the basis states according to the set 

{|1), |2), |3), . . ., |9)} = {\2),ai\ + l),al\ - l),aL| + 1), 
al\-l),aibl\0),albl\0),albl\0),albU0)}- (4) 

Initially the coefficients {Cq, C7, Cs, Cg} are prepared. Of these, the Hamiltonian couples in 
Case I Cq to C3 and Cg to C4 only; in Case II Cg to C3 only. In these subspaces the system 
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can be solved exactly, and performing a rotating wave approximation with respect to the 
frequency ujq + fli + fl2 we obtain in Case I 

Cg{t) = cos(Ait)C9(0) + i sin(Ait)C4(0) (5) 
Ce{t) = cos(Ait)C6(0) + I sin(Ait)C3(0); (6) 

in Case II only is valid. Choosing the interaction time such that Ait = vr, we find that the 
probabilities are restored in these subspaces. We are now left in Case I with a 5 dimensional 
and in Case II a 7 dimensional subspace to consider numerically. 

After the interaction, the state (|^) is available for measurements. In the ideal situation, 
the initial photons would have been restored to the radiation field. This is desired because 
the information resides in these photons, and they should be available for subsequent com- 
putational operations. We can ensure that they have been returned by observing that the 
atom is back in its ground state |0) by projecting the final state on this. After the interac- 
tion, the atom is available for inspection; a measurement on its state does no longer affect 
the outcome of the process. We write this state after an observation, |\l/o) = \0) {^l^^ out) , as 

+ ( -^e*''+-a+at + -^e*'^--a_aO (3^bUo). (7) 

a-p- J 

We have written the amplitudes and phases of the new coefficients as Cj^e*'^*^' G {— , +}). 

A measure of the efficiency of the process is the probability Pq = I (Ol^^out) A small 
value of Pq makes the process inefficient, but once the state |0) has been observed on the 
atom, the expressions in the brackets of (|^) give the effect on the state (a+a+ + a„al)|0) 
conditioned on the presence of the photons b^. on the lower transitions. These expressions 
contain the effect of the gating action of the system. In all cases investigated in this paper, 
however, Pq has been found to deviate from unity by less than 1%. The process is efficient 
as given. 

If the coefficients rjij = Cij/\aiPj\ in (|^) are close to unity, the interaction only adds the 
phases cpij] the polarization of the a"'"-field has been changed by the interaction. If we define 



the initial phases ip'^ = arg{a±) and ip\. = arg{P±), we denote the phase changes by 

^± = ((p+± + -ipl- ipt)/2 - ifl (8) 
Ayp± = (^+± - <^-± -^1 + <^^)/2. (9) 

We now write the final state (|^) in the form 

+e^^- e^^^- a+a^ + r/_^e-*^^-a„ al) (3^bl} |0). 

When we choose the initial coefficients a±,P± real, the phases (H)-® simplify; at the end 
of this paper we are going to discuss the inffuence of the phase on the gating performance. 

In Case I, the symmetry requires that if_^^ = Tp_ and A</)+ = — Ay9_ = Aip. In 

Case II, we assert that ~ ip ~ 0, which implies Tp_ ~ and A(p_ ^ 0. We 

may consider the 4-dimensional subspace {at^fot. |0), |0), at.6+|0), a+6+|0)}. Assuming 
now that all coefficients rjij are unity, we obtain in the symmetric case the ideal trans- 
formation Ui ~ e^'^Diag{e'^f,e-^^'f,e-^^'^,e^^^}. In the detuned Case II, we obtain 
Uii ~ Diag{l, 1, e^(^++^^+), e^(^+-^¥'+)}. This is a phase transformation of the bit denoted 
by a\. induced by the presence of the photon h\. 

We are now going to consider the performance qualities of the model system as a gated 
bit transformation. The input to the calculation is the initial state @). To begin we choose 
the "classical" case when only one of the input states is present. In the symmetric Case I, 
the choice of state is not important, c.f. [//, but for the Case II, we need to look at the 
states a^fo^^lO) and a^6^|0). First we choose to discuss the single input state al6^|0) with 
«_=/?+ = 1. 

As stated above, the interaction time is chosen such that t = vr/Ai; in the calculations we 
choose Ai = 1. For large detunings {lj2 oo) rj_^ approaches unity, but the phase shift A(y9 
goes to zero. In Case I, the numerical investigations show that we can retain rfi^ > 0.9 if we 
choose A2 > 5. For Ag = 5 we find A<y9+ ~ 10°. This is achieved with A2 = 1; larger phases 
can be achieved by increasing A2, but the restoring of the population suffers. For A2 < 1.5 we 
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can achieve A(/)+ > 15° and 77^^ > 0.75. The results can be illustrated in a graph plotting 
A(y9+ as a function of 7f_j^ with the detuning as a parameter. For the symmetric Case I, 
this is done in Fig.2a. As we can see, for A2 > 5, no dependence on detuning is seen. The 
corresponding results for Case II are shown in Fig.2b. Here the dependence on detuning is 
much stronger; however for large values of detuning, = 15 and = 30, we can reach 
A(p+ > 43° with _|_ > 0.9. Thus the operation of this gate is much more efficient as is to be 
expected. For larger values of A{^ the results tend to become independent of the detuning. 

We now choose to look at the case A2 = 2.5 and A2 = 30. For the Case I this gives 
A(y9+ ~ 10° and r]'^_^_ ~ 0.90. In Case II it gives A(/9+ ~ 10° and r]'^_^_ ~ 0.99. In order 
to see where the missing population goes in Case I, we plot the population of the states 
at.&+|0), at.|0), a+|0), a\^bL\0) and |2) in Fig. 3. At time 1 = 11, the population of aLb\_\0) is 
restored to 90% but the missing population is on the level a+fet. |0). This is mediated through 
the off-resonant transition | — 1) ^ |2) ^ | + 1) which proceeds at the effective Rabi rate 
(A2/A2) ~ 6.25/30. With time, this increases the population of the state + 1) as can be 
seen in Fig. 3; this increase is modulated at the rate Ai by the population pulsations on level 
at I — 1). This effect can be decreased by increasing A2 >> A2. In Case II, the population 
of the level aLb\_\0) is restored to better than 99% and the population on states a+| + 1) 
and albl\0) remain below 10 . 

After having described the "classical" inputs, where each 2-bit pure state has been treated 
separately, we now turn to consider the genuine quantum situation described by the input 
state d^). The performance of the system acting on this state is, of course, essential for its 
usefulness quantum computing device. 

An input consisting of a pair of two-level systems contains 4 degrees of freedom: the 
4 complex numbers involved loose two parameters to the over-all phase and two to the 
normalization conditions. It is still difficult to display the results of a 4 parameter input 
space, and hence we start by considering only real coefficients in Eq.(|^). The influence of 
the phases ^9^, y^^. will be discussed below. We are thus left with two real parameters, one 
for each input bit. We choose to display our results as functions of = 1 — for the two 
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cases 



lA) = ;^(i.UM_j|0) ^^^^ 
Ift) = (fit + ift*-) |0> 
We want to introduce a quality factor for the use of a system like this in computations. 
The performance is close to ideal, when the parameter rjij ~ 1. However, when either one of 
the input parameters ctj, l3j becomes close to zero, any minute value in the corresponding co- 
efficient Cij is likely to cause a large value rjij. Thus we want to consider the retention of that 
product aif3j which is the largest. A value close to unity here signals a good performance. 
To test this idea we consider the variables 

V-+ (a- > 0.5) ; ril_ (a^ < g.S). (11) 

Another measure of the efficiency of the process can be given by the retention of the ratio 
between the two components in Eq.(|^). This starts from and if retained the 
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should be close to unity. The retention parameter R for Case I and the inputs and 
1/32) is shown in Fig. 4a together with the corresponding quality factor in Eq.(0). In Fig. 4b 
the same parameters are shown for the asymmetric Case II. As we can see, the retention 
parameter R is at its worst about 70%; in Case II it is better than 90%. In Case I, the 
quality factor ([TI|) is good to within 90% and in the asymmetric Case II to better than 95%. 

Finally we want to return to the question of the influence of the initial phases. These do 
affect the outcome, but their influence seems to be smaller than the influence of the magni- 
tudes. We consider the achieved phase shifts as functions of the superposition coefficients a 
and (3. In Fig. 5 we plot the phase shifts A(p± against in the asymmetric Case II shown 
for and |/32). For we also consider the case when the initial phase (p^. is set to 
the value y?" = ir/A. The behaviour is close to ideal; in the range G (0.1,0.9), nearly 
ideal behaviour is observed, A<y9+ ^ 9.5° and |A<yO_| < 0.4°. The effect of the initial phase 



is small. In the symmetric Case I, the behaviour was found to be less optimal: we saw only 
a small difference for the two /5-states, but for o?_ in the range (0.1,0.9) the phase shift 
changed from 30° to 10°. Thus in Case I, the magnitude of the angle remains considerable 
but it does depend on the value of a. We have not carried out a systematic investigation 
of the influence of the phase factors; the results reported here indicate that they cause no 
drastic changes. If needed, their effects can easily be evaluated using the method presented 
here. 

As a conclusion, we discuss how well a quantum gate can be realized in our model. 
We choose to look at the Controlled-NOT gate, which changes the value of the target bit 
whenever the control bit has the value one. Based on the considerations above, we conclude 
that the asymmetric Case II is better suited to work as a gate. Its performance can easily be 
improved from the results presented above by increasing Ag^, and A2 in a suitable way. 
Here we use the parameters IS^^ — "^^i — 65, A2 = 6.85, Ai = 2, and t = n: this enables 
us to approximate the transformation Ujj to the accuracy 10^^ with a phase shift of 60°. 
This has to be applied three times in sequence in order to get a phase shift of 180°, which 
is needed for the Controlled-NOT gate. After performing suitable transformations between 
the circular and linear bases (see we obtain as the final result the Controlled-NOT 
transformation Cn'- 

'0.9956-^33° 0(10-3) 0(10-2) C(10-2) 
C(10-3) 0.995e-^33° 0(19-2) O{10-^) 
C(10-2) C(10-2) C(10-3) -0.997 
C(10-2) C(10-2) -0.997 C(10-3) 

The overall phases e-*^^" and —1 are irrelevant. We see that the Controlled-NOT gate can 
be realized in this case to the accuracy 10-^. 

The present scheme has been found to perform reasonably well as a computing device. 
It is naturally not good enough to be an element of a computer network of realistic size, 
but there seems to be no suggestion in the literature which satisfies this criterion. The 
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performance of our scheme can be improved by sequantial application of the b^- and a^- 
photons, with final restoration of the 6^-state by a third pulse. Such a scheme seems to 
require perfectly controlled pulses, which we regard as even more unrealistic than the model 
we have investigated. To implement our method in a multi-step computation we assume 
all initial information to be coded in a set of field modes residing uncoupled in the same 
cavity. During their coherence time, we shoot through the cavity volume a sequence of 
suitably chosen atoms which couple the photon pairs, i.e. perform the two-qubit operations. 
To affect all possible unitary transformations, the cavity has to be rather complicated, 
containing a suitable arrangement of A-plates to give access to all desired polarization states. 
Also the atoms have to be able to couple just the desired modes at each stage of the 
calculation. This and the restrictions imposed by loss rates and decoherence times pose 
extremely strict limitations on the computations possible. If several cavities are necessary, 
the dissipative effects on photons transferred between them raise further problems. However, 
such difficulties seem to afflict other schemes suggested too. Which one can be optimized 
the most remains an experimental challenge. 
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FIGURES 




FIG. 1. The 4-level system used in the gated Faraday effect 
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Figure 2a Figure 2b 




FIG. 2. The phase shift A(^+ as a function of rj^j^ for several values of the detuning A^, Ag^; 

some values of A2 used are marked. In all figures Fig.^a corresponds to Case I, Fig.T^^b to Case II 
where A(^ = ^2/2. 
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FIG. 5. The phase shifts A(/?± as functions of a?_, for (solid hnes) and |/32) (dotted hnes). 
The shift A(^_|_ is shown also for the case of a non-zero initial phase ^p'i^. (Case II) 
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